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We study the properties of fermion correlators in a boundary theory dual to the Reissner- 
Nordstrom AdSd+i background in the presence of a bulk dipole (Pauli) interaction term with 
strength p. We show that by simply changing the value of the parameter p we can tune con- 
tinuously from a Fermi liquid (small p), to a marginal Fermi liquid behavior at a critical value of 
p, to a generic non-Fermi liquid at intermediate values of p, and finally to a Mott insulator at large 
values of the bulk Pauli coupling. As all of these phases are seen in the cuprate phase diagram, the 
holographic model we study has the key elements of the strong coupling physics typified by Mott 
systems. In addition, we extend our analysis to finite temperature and show that the Mott gap 
closes. Of particular interest is that it closes when the ratio of the gap to the critical temperature 
is of the order of ten. This behavior is very much similar to that observed in the classic Mott 
insulator VO2. We then analyze the non-analyticities of the boundary theory fermion correlators 
for generic values of frequency and momentum by calculating the quasi-normal modes of the bulk 
fcrmions. Not surprisingly, we find no evidence for the dipole interaction inducing an instability in 
the boundary theory. Finally, we briefly consider the introduction of superconducting condensates, 
and find that in that case, the fermion gap is driven by scalar-fermion couplings rather than by the 
Pauli coupling. 



I. INTRODUCTION 

Holography can offer unprecedented insight into the 
dynamics of strongly coupled systems. In the recent past, 
it has become clear that the domain of applicability of 
holography goes beyond high energy physics and includes 
strongly-correlated sytems in condensed matter physics, 
as well (see [1-4] for reviews). Since in most cases of in- 
terest, we do not possess a microscopic understanding of 
the held theory dynamics under study, a phenomenologi- 
cal point of view (the so-called "bottom-up" approach) is 
taken where a minimal gravitational setup is devised for 
analyzing a specific strong-coupling feature of a system. 
In some situations where string or M-theory completion 
of a bottom-up construction is known, one may wonder 
how a bottom-up result is modified in a top-down ap- 
proach. For example, a quantity of interest in strongly- 
correlated condensed matter systems, which can easily be 
computed using holography, is the fermion spectral func- 
tion (which is proportional to the imaginary part of the 
fermion retarded two-point function) . There has recently 
been much discussion about this quantity in the holo- 
graphic literature [5-9] where analyzing the Dirac equa- 
tion for a charged probe fermion propagating in a grav- 
itational background (usually a charged black hole), one 
can show that the retarded two-point function of the dual 
fermionic operator in the boundary theory shows a vari- 
ety of unexpected emergent phenomena. In top-down ap- 
proaches to holographic systems, the fermions are gener- 
ically coupled to gravity and gauge fields in a variety of 
ways, beyond minimal coupling. It is certainly desirable 
to analyze how such non-minimal bulk couplings modify 
the fermion spectral function in the boundary theory, or 
may lead to new interesting emergent phenomena. 

Recently, it was shown [10-f 3] that there are consistent 
truncations of ten and eleven dimensional supergravities 



to five and four dimensional bulk theories that possess an 
interesting class of gauge interactions and charged mat- 
ter, allowing for novel condensed matter physics phenom- 
ena, such as superconductivity, to be explored in a consis- 
tent top-down approach [14-16]. The fermionic sector of 
these truncations has also been worked out in [17, f8] 
where a number of generic couplings for the fermions 
(with possible applications to strongly-correlated con- 
densed matter systems) have been realized. 

Motivated by these studies, we considered in [f9] a 
generic non-minimal fermion coupling in which a spin- 
1/2 fermion couples to the gauge field through a dipole 
(Pauli) interaction of the form F a t,'ipT ab ip. In fact, we just 
considered the simplest possible setup, in which a spin- 
1/2 fermion propagates in the background of a Reissner- 
Nordstrom AdS^+i black hole. We found that as one 
changes the strength of this interaction, spectral weight 
of the dual fermionic operator is transferred between 
bands, and beyond a critical value of the dipole coupling, 
a gap emerges in the fermion density of states. We then 
concluded that a possible interpretation of this interac- 
tion is that it drives the dynamical formation of a (Mott) 
gap, in the absence of continuous symmetry breaking. In 
the Hubbard model, the Mott gap forms in d > 1 once the 
on-site interaction U exceeds a critical value in the half- 
filled system. Upon doping, spectral weight shifts from 
high to low energies. Consequently, we argued that the 
strength of the dipole interaction mimics the combined 
effects of doping and the on-site interaction strength U. 

In this paper, we continue our study of the dipole in- 
teraction in more detail. We investigate the existence of 
Fermi surfaces as the dipole coupling p is varied. For 
the range of parameters considered in this paper, we find 
that there is no Fermi surface above a certain value of p. 
This is the range of parameters considered in [19], where 
the boundary theory exhibits two main features of Mott 
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insulators, a dynamically generated gap (in the absence 
of continuous symmetry breaking) and spectral weight 
transfer. In addition, for the values of the dipole cou- 
pling p for which there is a Fermi surface in the boundary 
theory, we find that at a critical value of this coupling, 
marginal-Fermi liquid behavior ensues. Consequently, we 
are able with our model (see Figure 1) to describe at one 
extreme, the Mott insulating state with a dynamically 
generated gap, a transition to a marginal Fermi liquid 
(MFL) and at the other, a Fermi-liquid-like regime in 
which the excitation spectrum scales linearly with the 
frequency. As all of these regimes are accessed simply 
by changing the value of the dipole coupling p, this sug- 
gests a direct parallel between p and the hole-doping level 
in the high-temperature copper-oxide superconductors 
(hereafter cuprates). In the cuprates, the strong electron 
correlation physics ends at a value of doping (typically 
optimal doping) where marginal-Fermi liquid [21] behav- 
ior ensues. Perhaps the occurence of the MFL state of 
the probe fermions in our holographic setup is an indi- 
cation that this phenomenological model is ultimately a 
robust feature of the transition from strong to weakly 
interacting physics in doped Mott systems. 

Having discussed the pole structure of the boundary 
theory fermion (retarded) correlators for non-zero values 
of p at small frequency, we analyze the non-analyticities 
of those correlators for generic values of frequency and 
momentum. We do this by numerically calculating the 
quasi-normal modes of the bulk fermion in the (extremal) 
Reissner-Nordstrom AdS^+i background. Following the 
motion of the poles in the complex frequency plane as 
a function of momentum, we compute their dispersion 
relations and, for the range of parameters considered, 
confirm that all of the poles stay in the lower half of the 
complex frequency plane, for all momenta. Hence, as 
expected, turning on a non-zero bulk dipole coupling in 
our set up does not cause an instability in the boundary 
theory. 

We consider our setup at finite temperature and find 
that as the temperature increases the gap closes and, 
moreover, the critical temperature for which this happens 
is much less than the value of the gap. In this regard, the 
closing of the gap induced by temperature parallels what 
one obtains in the classic Mott system VO2 [20] , namely, 
the gap closes at a temperature much lower than the gap. 

The paper is organized as follows. In section II we con- 
sider a bulk fermion in the Reissner-Nordstrom AdS^+i 
black hole background and couple it to the U(l) gauge 
field non-minimally through a dipole interaction with 
strength p. We then derive the Dirac equations and 
rewrite them as flow equations which will be more con- 
venient for numerically calculating the boundary theory 
fermion correlators. In section III we investigate the ex- 
istence of Fermi surfaces as a function of p by solving the 
Dirac equations at zero frequency. In section IV, we first 
discuss the small-frequency behavior of the poles of the 
boundary theory fermion (retarded) correlators when p is 
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FIG. 1: A cartoon of the phase diagram of the boundary theory 
considered here. MI indicates a Mott insulator, a phase with a 
gap in the absence of symmetry breaking. NFL denotes non-fermi 
liquid behavior which is distinct from the gapped spectrum of a 
Mott insulator. MFL indicates marginal-Fermi liquid behavior in 
which the "electron" self energy scales as uj logo; at T = 0, and a FL 
(Fermi liquid) regime in which the dispersion is linear in frequency. 
The tuning parameter in this model is the Pauli (dipole) coupling. 
Similar behavior is obtained in the non-superconducting features 
of the cuprate materials by tuning the hole-doping level x. 

non-zero. We then analyze the non-analyticities of those 
correlators for generic values of frequency and momenta. 
In section V we study the effects of temperature in our 
holographic setup. Finally, in section VI, we discuss the 
relevance of our work to cuprate phase diagram and con- 
clude with open questions as well as extensions for future 
work. In particular, we contemplate the extension of our 
results to superconducting backgrounds, which also do 
not possess a finite ground state degeneracy at zero tem- 
perature. 



II. BULK ANALYSIS 

As we alluded to above, we consider just one form of 
non- minimal coupling, in which a spin- 1/2 fermion is cou- 
pled to the gauge field through a dipole interaction 1 of 
the form F a t'tpr at 'ip , and propagates in the background 
of a Reissner-Nordstrom AdS^+i black-hole (hereafter, 
denoted by RN-AdS^+i). Thus, we consider the bulk 
Lagrangian 

V '—gi't/j(Jp — m — ipf)"4>, (1) 
in d + 1 > 4 dimensions where 

U) = efr c (d M + \u$T ah - iqAjA , 
f=\v ah e™e%F MN , ( 2 ) 



In even bulk dimensions, there is a similar interaction which 
includes a F$. We will not consider this interaction in this paper. 
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with e^ 1 and uj^f being the (inverse) vielbein and the 
spin connection, respectively. We denote the bulk co- 
ordinate indices by capital letters M, N, ■ ■ ■ = {t, x l , r} 
while the tangent space indices are denoted by a, b, ■ ■ ■ = 
{t,x l ,r}. We will reserve the Greek indices fi,v,... 
to denote boundary coordinate directions. We use 
Dirac matrices r*, ...,T- satisfying the Clifford alge- 
bra {r a ,r fc } = 2rj ab . Also, T ab = |[r a ,r 6 ]. In what 
follows, we will rescale p — > pL/(d — 2) for convenience. 

The RN-AdSd+i background has a metric and a gauge 
connection which can be written 



ds 2 



L 2 



[-.f(r)dt 2 + dx 2 ] + 



L 2 dr 2 

r 2 J( r ) ■ 



A = pt 



1 



r 



dt, 



(3) 
(4) 



where 

m = 1-^(7) 

/ d- 

H = 



)2 ^ro^(d-i) 
r 



d - 1 \ 1/2 Qr 



2d -4 



L 2 



M = l + Q 2 , (5) 



with r being the horizon, given by the largest real root 
°f f( r o) = 0. The temperature T of this (black hole) 
background is given by 



T= (d-2) 



d 



4ttL 2 V d - 2 



Q' 



(6) 



From the above equation, one notes that the RN-AdS^+i 
black hole is extremal when Q 2 = d/(d — 2) while the 
density and entropy remain finite. Since the background 
is invariant under A t — > —At, without loss of generality, 
we can choose fi, or equiv alcntly Q, to be positive. Thus, 
we can take < Q < yjd/{d— 2), where the equality 
corresponds to extremality. 

To analyze the Dirac equations of the bulk fermion, 
we find it more convenient to go to momentum space 
by Fourier transforming ip(r, x^) ~ e tk ' x ip(r, k^), where 
fc M = (00, k). The Fourier transform of the Dirac operator 
ifi is of the form 



o f(r) d 
4/(r) 2r 



L 



— i- 



VW) 



while 



id - 



f=(d-2)^(^\ 
r \ r J 



qMr)} 



i — k ■ f , 
r 



d-l 



(7) 



(8) 



To decouple the Dirac equations, we introduce pro- 
jectors r± = \{l ± Y^Vi-k ■ f) and write ip±(r) = 
r d / 2 f{r) 1 / i T±ij)(r). Without loss of generality, we set 



k\ = k and k^\ = 0, and take the basis 

r , , <T 3 ®1 \ t _ /t<7i<g)l 

" ! -cr 3 ® 1 i ' I io-i ® iy ' 



1 _ f-a 2 ®i 

0-2 ® 1 y ■ 



(9) 



where cr/s are the Pauli matrices, and 1 is a 2 2 _ 

d — 4 

dimensional identity matrix for odd values of d, and 2~- 
dimensional for d even. Note that by choosing k\ = k and 
kijti — the rest of the gamma matrices do not appear 
in the Dirac equations. So, we did not bother to include 
those in (9). One then finds 



^c^2 



0J + fiq\ 1 



„d-2 



- °\ W—rji ± k \ip± - a 3T mip±. (10) 



We see that the Pauli coupling modifies the appearance 
of k in the above Dirac equations. To see the effects of 
p more clearly, consider the solutions of the Dirac equa- 
tions (10) in the asymptotic and near horizon regimes. 
Asymptotically, the solutions behave as 



<P±(r,LJ,k) = a±(LO,k) r mL (J 



+ b ± {oj,k)r- mL (J 



1 + - 
1 + 



(11) 



The effect of p asymptotically is to modify the subleading 
terms. For me [0, \) both terms in (11) are normalizable 
and one can choose either a± or b± to be the sources for 
the dual fermion operator in the boundary theory. In this 
paper, we take m E [0, |) and consider the conventional 
quantization where a± are the sources. Thus, the dual 
fermion operator has dimension A = | + m. Choosing 
in-falling boundary conditions near the horizon results in 
a retarded correlator of the form 



G R (u>,k) 



G + {u,k)l 
G_ 





(u,k)t 



(12) 



with G±{oj,k) = b±(uj,k)/a±(uj,k). Note that the Dirac 
equations (10) imply G + (uj,k) — G_(w,—k). 

When the background is extremal, f(r) has a double 
zero at the horizon, f(r) ~ d(d — 1)(1 — r /r) 2 + ■ ■ ■ , 
and this fact makes taking the limit of ui — > of the 
equations (10) near the horizon subtle. To take care of 
the subtlety, one realizes [8, 22] that near the horizon (in 
which the geometry approaches AdS2 x R d_1 for T = 
0) the equations for ip± in (10) organize themselves as 
functions of £ = LoL 2 /(r — r ) with L 2 = Lj \fd{d — 1) 
being the radius of AdS2 . The coordinate ( is the suitable 
radial coordinate for the AdS2 part of the near horizon 
region, and in this region, we can write ip± in terms of £ 
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and expand in powers of w as follows 

v>/± (0 = v4°2 (0 + w vfi (0 + w vfi (0 + 



(13) 



Now, substituting (13) into (10), we find that to leading 
order 



hi 

c 



TflC 3 



P , kL\ 
c d -± 



r o J 



4°2(0 



- i (72 ( 1 + 



qed 

c 



V4°i(0, 



(14) 



where e d = l/ v / 2d(d- 1), and c d = 1/ [(2d - 4)e d ]. 
Equations (14) are identical to the equations of motion 
for massive spinor fields [8] with masses (m, m+) and 
(m,m_) in AdS2, where m± are time-reversal violating 
mass terms, with the identification 



m± 



p kL 

Cd-r ± • 

L r 



(15) 



Thus, i/)j°2(Q are dual to spinor operators 0± in the IR 
CFT with conformal dimensions 5± = + | where 



' rn 



m\±= m 2 



Cd T ± — 
L r J 



(16) 



We see that turning on p modifies the scaling in the in- 
frared in an important way - effectively, the momentum 
is pushed up and down by p. We will explore the details 
of this in what follows. 

One can write a formal expression for the fermion re- 
tarded correlator (12) at low frequency in terms of the 
retarded Green functions of the IR CFT spinor opera- 
tors 0±. As shown in [8], such a formal expression is 
extremely useful in analyzing the small ui behavior of 
the boundary theory Green functions. This is done by 
matching the inner AdS2 and outer AdS4 solutions in 
the so-called "matching region" where the £ — > and 
uj/( — >• limits are taken. In so doing, one finds that the 
coefficients a±(co,k) and b±(ui,k) in (12) are given by 



a±(iv, k) = 


a± 




+ 


~(0) 

a± 


+ u~a2 ] 


b±(w,k) = 






+ 




+ ub¥ 


where a±\ a±\ b±^ 


and 



with 



c±{k) 



T(-2v ± )T(l + v ± -iqe d ) 
T{2v ± )T{l- V± -iqe d ) ' 
(m + im±)L 2 — iqe d — v± 
(m + im±)L 2 - iqe d + u± ' 



(20) 



Note that the expressions (17) are (18) are not valid when 
2is± is an integer. In such cases there would be additional 
terms like u n \oguj (with n being a positive integer) on 
the right hand sides of (17) and (18). 

In order to obtain G± (w, k) for generic values of ui and 
k, one must solve the Dirac equations (10) numerically. 
For numerical purposes, it is convenient to work with 
dimcnsionless quantities. So, we rescale r, lo and fc in the 
Dirac equations (10) by defining 



r u, 



k^^k. 



(21) 



It is also more convenient to convert the Dirac equations 
(10) into the so-called flow equations [6, 23]. For that, 
we first write V± = {fl± , a ± ) and define £± = (i± /a± , in 
terms of which the Dirac equations (10) then reduce to 
the non-linear flow equations 

u 2 ^J(u)d u U = -2(mL)uU 

+ [v_(u)T A] + [«+(«) ± *] (22) 



where 

v±(u) 



[u + Qq(l -u 2 - d )] ±Qpu 2 - d . (23) 



To obtain the retarded Green functions of the boundary 
theory operators, one has to choose the infalling bound- 
ary condition at the horizon for the (dual) bulk fields 
[24, 25] . Expressed in terms of £± , the infalling bound- 
ary condition for ip± = (f3± ,a±) at the horizon translates 
into 



1 (24) 
(mL 2 - v±) / (qe d + fh±L 2 ) u> = 0. v ' 



e±(« = i) = 

The matrix of Green functions (12) then takes the form 



O (u 2 )] gftu,), (17) Gr ( U , k) = lim e- 2mL {^f ^ 



(25) 



(18) 



are all functions of k and 
can, in principle, be determined numerically. Also, in 
the above expressions, Q^{uj) denote the retarded Green 
functions of the dual IR CFT operators 0± which are 
given by [8] 



(19) 



where one picks the finite terms as e — > 0. Up to normal- 
ization, the fermion spectral function is defined by 

A(u},k) = TrTmG R (u>,k). (26) 



III. CONTINUUM AND BOUND STATES 

In this section, we will study the effects of small and 
negative values of p. As we will see, in this regime, 
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there is a (non-Fermi-liquid) Fermi peak whose proper- 
ties change as we vary p. In this regime, there is some 
similarity to the properties of Fermi surfaces studied in 
[6, 8]. As we change parameters, the scaling dimensions 
change, and we can pass from a non-Fermi liquid (NFL) 
to a marginal Fermi liquid (MFL) and on to fermions 
which have some resemblance to Landau Fermi liquids 
(FL). However, we will see clearly that there is a positive 
value of p beyond which the Fermi peak ceases to exist. 
In what follows, we will mostly set m = and 5 = 1 (as 
we vary p), although similar results hold for a range of 
these parameters. 

We begin by focussing on the regime that has been 
called log-oscillatory in [6, 8], in which Fermi peaks do 
not occur. There are some important changes whenp ^ 
that we will explain below. When q 2 > 2m 2 L 2 , there ex- 
ists a range of momenta k G X + for which the dimension 
of the IR CFT operator 0+ becomes imaginary. Sim- 
ilarly, for k G I- the dimension of 0_ becomes imagi- 
nary. Here, we have defined Z± = (T c d,p—k ,TCdP + k ) 
with k = yj (q 2 /2) — m 2 L 2 . (Figure 2 shows plots of 
I± versus p.) Consequently, ImG±(0, k) is generically 
non- vanishing for k G T± 1 respectively; see Figure 3 for 
plots of ImG±(0, k) as a function of p. At p — 0, one 
has I + — I_ = I. This case was analyzed in [6, 8] 
where it was found that for fcel and for small u>, both 
ImG±(cLi, k) are periodic in logo;, with the same period. 

At p = 0, the range of momenta for which Im G± (lo, k) 
become log-oscillatory at small oj is the same for each, 
namely k G I. This degeneracy does not persist 
for non-zero p, hence the fermion spectral function 
A(lu, k) will also have non-oscillatory components. For 
p G [— k a /cd, k /cd] — {0}, both lmG±(uj,k) show log- 
oscillatory behavior (with different periods, though) only 
for k el + fll_. For \p\ > k /cd where l + nl_ = 0, one 
finds that in the regime where ImG_(oj,/c) shows log- 
oscillatory behavior, ImG+(w,/c) is not oscillatory and 
vice versa. 

For real v k , the boundary conditions for £±(u = 1) 
at co = are real. Since the equations (22) are real, 
one deduces that ImG±(0, k) — 0. Thus, ReG±(0, k) — 



(0) /ai 0) 



G±(0,k) = b 

would be given generically by the zeros of a±' . Each zero 



ReG±(0, k) may have poles which 

(o) 



of cl± defines a Fermi momentum kp, given that b K ±' do 
not vanish as k — > kp. Since G_(w, k) — G + (u;, — fc), van- 
,(°) 



(0) 



ishing of a}" ' at some k = kp implies that a/? vanishes 
at k = —kp. So, in order to find fcp, we can just analyze 

the zeros of . From the asymptotic behavior of tp- 
which is given in (11), together with (17) and the defini- 
tion of ifL = (/?- , a 
a-[u,k) = a { " } u mL 



), it is easy to see that, at u> = 0, 

• • • as u — y oo. (Also, note that at 

u) = 0, P-(u,k) =6 (0) M _mL + --- asw^oo.) Sothe/c F 's 
define a set of momenta for which, at lj = 0, tp- (u, k) be- 
comes normalizable (a 'bound state') as u — > oo. To find 
the fcp's, we analyze the equation for a_(u,fc) as follows. 
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FIG. 2: Plots of Z± versus p (for d = 3, = 3/2 and q = 1). The 
red band depicting I_ is where Im G— (uj, k) becomes oscillatory at 
small u). The blue band (1+) shows the region where ImG+(u, k) 
is oscillatory (at small u)). 



Plugging = (j3- , a_) into the Dirac equations (10), 
one obtains a set of two coupled linear differential equa- 
tions for a_ and (3- . Setting to = and decoupling these 
two equations, we obtain 



v°_(u) + k 

u 2 VK 



u) I u 

On 



v+(u) — k 

2 vm 



v+(u) — k u \ v°_(u) 



d u a- 



(27) 
(28) 



where we have set m = for convenience. In (27), the su- 
perscript "0" on v± (u) indicates that we have set lj = 
in (23). Once again, we set d = 3 and q = 1 in what 
follows. As we vary p, we look (numerically) for the mo- 
menta kp for which a_(it = oo,kp) — 0, given an ap- 
propriate boundary condition for a_(u,u; = 0,k) at the 
horizon. Indeed, solving (27) near the horizon, one easily 
obtains that a_(u,k) ~ f(u) ±Uk ^ 2 as u 1. Because 
by assumption we are in a regime where v7 is real and 

positive 2 (and, in fact, generically irrational), /(u) - '^ I 2 
blows up as ii -> 1. Thus, a_(u l,k) ~ f(u)" k I 2 is 
the regular horizon boundary condition that should be 
chosen. 

In Figure 4 we have plotted such values of hp as a 
function of p. Starting with negative values of p (while 
keeping q = 1 fixed) , kp increases as we raise p causing it 
to move towards the boundary of the oscillatory region 
I_ . As p approaches l/\/6 from below, kp approaches \/2 
(in units of r /L 2 , from below). The blue dots in Figure 



2 For vZ = 0, one finds that ct-(u — > l,fc) = a(l-| ) + 61og(u — 

1)(1 + - ■ ■ ) where the dots represent terms which vanish as u — > 1, 
and a and b are some constants. In order for o_ (u, k) not to blow 
up at the horizon, one should then choose b = 0. 



6 



1.5 

1.0 

+i 

1 0.5- 
0.0- 



(a) 



1.5 

1.0- 
O 

1 0.5- 
0.0- 



(b) 



1.5 
o 1.0- 
J 0.5- 



(c) 




FIG. 3: Plots of Im G_ (0, fc) (solid line) and ImG+(0, fc) (dashed line) for (a) p = 0, (b) p = 0.1 and (c) p = 1.8. We set d = 3, = 3/2 
and g = l. Similar plots can be obtained for negative values of p by switching the solid lines with the dashed lines. Focusing on positive 
p, we see that the maximum value of ImG±(0, fc) increases as p — > after which (namely, for p > l/\/6) it rapidly decreases. 



4 show the location of kp's versus p and the red band 
depicts the oscillatory region. There is a single Fermi 
surface for each p as we increase p up to p = 1/ v6. We 
have explicitly checked that j3— (u — co, fc — kp) does 
not vanish, so k — k-p are genuine poles of ReG_(0,fc). 
At p = k F = and as a result v k ^ vanishes 

(recall that d = 3, m = and q = 1). (At this point, 
and in fact at any point in which 2f fcp G Z, the analysis 
should be more carefully done, as logarithms must be 
included.) For p > l/v6, we do not see a Fermi surface 
as a_(it = co, fc) does not vanish outside the oscillatory 
region. We have checked this numerically up to p = 10, 
and, given the observed behavior of a_ (u = co, fc), we do 
not expect it to change as we increase p further. Indeed, 
Figure 5 shows plots of a—(u — co, fc) versus fc for sample 
values of p. We have also plotted /3-(u — oo, fc), shown 
by the red curves in Figure 5. 

Following [8], the excitations around these Fermi sur- 
faces can be analyzed. Using (17) and (18), near k = kp 
and at small ui, G_(w, fc) takes the form 



G-(u,k) 



G) fe aL 0) (fc F )fc± + wa w (fc F ) + a^>(k F )g- (w) ' 

(29) 



& W (M 



where fcj^ = k—kp. Suppose the denominator in (29) van- 
ishes at some w*(fc) = Rew*(fc) — zlmcj^fc). For p = 0, 
the dispersion relation, the width and the residue of the 
pole were worked out in detail in [8] . Parts of the data in 
these quantities (such as the scaling of the dispersion re- 
lation) come from the IR CFT (or, equivalently, the AdS2 
part of the near horizon geometry) , and other parts (such 
as Fermi velocity) from the UV physics. For non-zero p, 
the data which come from the IR CFT will be slightly 
modified according to (16) whereas the data coming from 
the UV physics could be substantially modified. 

For -0.53 < p < l/v/6, wc find that 1/2 > > 0. As 
a result, the small u) excitations around fc = kp will have 
a (non- Fermi liquid) dispersion relation Rew*(fc) oc fc^ 
and a width Imw»(fe) oc fcr[ where z = \/{2v7). Thus, 
for this range of p, Ima;*(fc)/Rew*(fc) does not vanish as 
lj — > 0, implying that these excitations are not stable. 
Note that the residue at the pole is given by Z oc k^ 1 . 




FIG. 4: kp's (shown by blue dots) versus p. For p < l/v / 6 there 
is a single Fermi surface for each p. For p > 1/ \/6 we do not find 
Fermi surfaces. The orange band shows the oscillatory region I_ . 



At p — —0.53, =1/2 and the excitations near the 
Fermi surface are of the marginal Fermi liquid type. For 
-1.54 < p < -0.53, 1 > > 1/2, hence, the small w 
excitations around fc = fcp will have a linear dispersion 

relation Rew*(fc) oc k± and a width Imw*(fc) oc kj_ F . 
So, for this range of p, these excitations are stable as 
Imoj*(fc)/Rea;*(fc) — > as u> —> 0. It is in this sense that 
we refer to this region as the Fermi liquid. Also, we found 
that kp goes through zero at p = —1.317, signifying that 
the excitations change over from 'particle-like' to 'hole- 
like'. Consequently, we find that simply by varying p, we 
can tune from a Fermi liquid, p < —0.53 to a marginal 
Fermi liquid at p = —0.53 to a generic non-Fermi liquid 
for —0.53 < p < and finally to a Mott insulator 

for p sufficiently large. Precisely how the system behaves 
for large values of p will now be addressed in the quasi- 
normal mode analysis. 



IV. QUASI-NORMAL MODES AND STABILITY 

In this section we analyze the poles of G±(lu, fc) and, 
in particular, discuss how they move in the complex u>- 
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FIG. 5: Plots of a^(u = oo, k) (black curves) and 0—(u = oo, k) (red curves) versus k for (a) p = —0.4, (b) p = 0, (c) p = 0.2, (d) p = 1. 
The orange strip in each plot shows the oscillatory region I—. The plots are generated for d = 3, m = and 9=1. By k —> —k, similar 
plots could be obtained for a+(u = oo, k) and = oo, k). 



plane as we vary k. Since G+(w, — k) = G_(w,fc), with- 
out loss of generality, we can just focus on the poles of 
G-{ui, k). We denote the poles of G-(uj, k) by uj*(k). At 
small uj, the poles of G_(w,fc) can be worked out semi- 
analytically. Indeed, for p — 0, the small uj poles were 
worked out in detail in [8] for k £ 1 as well as k close to 
fcp where it was argued that such poles are all located in 
the lower half of the complex w-plane. Since the argu- 
ments of [8] are mainly based on the data coming from 
the IR region (the near horizon AdS2 region of the back- 
ground), they can easily be extended to non-zero values 
of p, where results similar to those in the case of p = 
are obtained. For example, at small uj and for k G I—, 
where ImG_(w, k) is oscillatory, the poles (for a fixed k) 
are exponentially separated on a straight line which is in 
the lower half cj-plane. (The line is diagonally oriented, 
i.e. it ends on the uj = branch point.) Equivalently, 
for those values of p for which there exists a Fermi sur- 
face, the small uj poles near k — fcp are all located in 
the lower half w-plane. In particular, if < 1/2, then 
uj*(k — > fcp) as a function of k follows a straight line in 
the lower half of the complex w-plane. The semi-analytic 
arguments of [8] , and their generalizations to non-zero p, 
are applicable only for small w, and for k £ I±, or when 
k is near ftp. Nevertheless, on general grounds, one ex- 
pects the poles to be located in the lower half w-plane 
beyond the small uj regime (and, of course, for all values 
of k). To find the poles beyond the small uj regime, one 
is usually forced to do numerics which, in the context of 
the AdS/CFT correspondence, involves performing some 
quasi-normal mode analyses in the bulk. 



Besides isolated poles, G±(uj,k) at zero temperature 
will have a branch cut (at least for small uj) which 
could be understood without doing the numerics. Note 
that since ~ u 2Uk appears in the expressions for 

a± (uj, k) and b± (uj, k) in (17) and (18), G±(uj, k) will have 
a branch point at u> = for generic values of k (where 
2vt are irrational), and a branch cut, which we take to 
be extended in the negative imaginary axis. For those 
values of k for which 2v^ g Z, there is still a branch cut 
which is due to the appearance of logarithmic terms of 
the form w"logw (n £ Z) in the expressions for a±(uj, k) 
and b±(u, k). The branch cut seems to be a distinctive 
feature of the two-point retarded correlators of operators 
in the zero temperature d-dimensional boundary theory 
dual to the extremal RN-AdS^+i background. Indeed, 
the branch cut was observed explicitly in the correlators 
of scalar and spinor operators in [8, 27] as well as the 
conserved currents in the shear and sound channels in 
[28, 29]. As we will see below, this branch cut appears in 
our quasi-normal mode analysis. At finite temperature, 
however, the branch cut dissolves into a series of isolated 
poles on the negative imaginary axis. 

Generically, Gn(uj,k) will have poles whenever 
a±(uj,k) = 0. In the context of the AdS/CFT corre- 
spondence, this problem could be addressed by comput- 
ing the quasi-normal modes of "0± m the RN-AdS^+i 
background, which are solutions to the Dirac equations 
(10) subject to the boundary conditions that they are 
infalling at the horizon and normalizable asymptotically. 
Except in very special cases, the generic values of the 
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quasi-normal frequencies are usually computed numeri- 
cally. We use the so-called Leaver's method [26] for this 
purpose. For concreteness, we take the boundary the- 
ory to be (2+l)-dimensional, i.e. d = 3. Our analysis 
can straightforwardly be extended to larger values of d. 
Also, as in the previous discussions, we consider m = 
and (7 = 1. 

Substituting tp± = (fi±,a±) in the Dirac equations 
(10), and setting m = 0, one finds 

u 2 d u (3± = [v- (u) Tk]a±, (30) 

u 2 ^fj(rf>d u a± = -[v+(u)±k]0 ± . (31) 

The equations for a± are the relevant equations for ob- 
taining the quasi- normal frequencies of tp± . Squaring the 
above equations, the decoupled equations for a± are eas- 
ily obtained 



u_(u)=Ffc \v + (u)±k 



d u \a± = -a±. 



(32) 



As we alluded to above, without loss of generality, we 
can focus on the quasi-normal frequencies of -0- and just 
analyze the equation for a_ in (32). 

In what follows, we switch to a new radial coordinate 
z = 1/u which is more convenient for doing the numerics 
in this section. In terms of the new radial coordinate, 
the horizon is at z = 1 and the asymptotic boundary at 
z = 0. The equation for a_ in (32) then becomes 



d 



(z) + k \v+(z) 



(33) 
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FIG. 6: (a) ImG-(w, k) as a function of uj for k = 2. (b) The 
quasi-normal frequencies of a_ for k = 2. d = 3, p = 5, q = 1 and 
m. = in both plots. Also, M = 250. 



where ^4~ p (w, k) are the elements of a (M+l) by (M + l) 
matrix A~(u>,k). The quasi-normal frequencies (for a 
fixed k) are then the solutions to 



det A~(u),k) = 0. 



(37) 



To compute the quasi-normal modes of ip-, the be- 
havior of a_ should be infalling at the horizon and nor- 
malizable at the boundary. As mentioned above, we use 
Leaver's method [26] to compute the quasi-normal fre- 
quencies. For that, we first pull out the leading behavior 
of a_ at the horizon as well as the boundary and write 

a_(z)=e'^)/(z) _ '( f+ ^)zL(z). (34) 

Note that ct-(z = 1) is a constant which could be set 
equal to unity as the equation for a_ is homogeneous. 
Next, we write ct_(z) as a power series in z around a 
point zo = 1/2 (so that the radius of convergence of the 
series covers both the horizon and the boundary) 

M / 1 \ m 

<M Z ) = "™( W ' fc ) [ Z ~ 9 J • (3 5 ) 

Substituting (35) and (34) into (33), one obtains 

M 

Y / A mp (o J ,k)a m (oj,k)=Q, (36) 

m=0 



The bottom plot in Figure 6 shows the quasi-normal 
frequencies of a_ for p = 5 and k = 2. To generate this 
plot we set M = 250. Due to space limitations, the plot 
only shows a handful of the quasi-normal frequencies. 
As M is increased, the poles located along the negative 
imaginary axis become closer to one another, suggesting 
that their existence is due to taking M to be finite, and in 
the limit of M — > oo they should indeed form the branch 
cut we mentioned earlier. On the other hand, increasing 
M does not seem to change the qualitative behavior of 
the poles which are oriented almost diagonally on each 
side of the negative imaginary axis. Notice that the poles 
are all located on the lower half w-plane. As may be seen 
from the plot, the branch cut bends to the right for large 
negative values of Imw. This behavior is different from 
the cases studied in [28, 29] where the unbroken parity 
symmetry of the boundary theory forces the branch cut of 
the retarded correlators to stay on the negative imaginary 
axis. Indeed, the bending of the branch cut is similar to 
what was observed in [27] for the retarded correlators of 
charged scalar operators in the presence of a magnetic 
field. The top plot in Figure 6 shows ImG_(cj, k = 2) as 
a function of u. The location of the peak on the left hand 
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FIG. 7: The top two plots show the dependence on k (dispersion relation) of the real and imaginary parts of the first five quasi-normal 
modes (depicted in Figure 6(b)) on the left hand side of the negative imaginary axis. The bottom two plots show the dispersion relation of 
the real and imaginary parts of the first five quasi-normal modes of Figure 6(b) which are on the right hand side of the negative imaginary 
axis. The plots are generated for d = 3, p = 5, m = and q = 1 and M = 250. The red data corresponds to the mode closest to the real 
axis in the complex aj-plane. 



side and the bump on the right hand side match quite well 
with the two quasi-normal frequencies which are closest 
to the real axis. It is apparent that all of the other quasi- 
normal modes are relatively wide and individually have 
small residue. 

The dispersion relation, cj*(fc), of the quasi-normal 
frequencies shown in Figure 6(b) can be computed nu- 
merically by following their motion in the complex oj- 
plane as k is varied. As shown in Figure 7, the two 
quasi-normal frequencies which are closest to the real 
axis have markedly different dispersion relations than the 
rest (higher resonances). The effects of these higher res- 
onances become important as one goes beyond the small 
frequency and momentum approximation. For large 
the imaginary part of the quasi-normal frequencies is ap- 
proximately constant (and negligible compared to the 
real part). In contrast, for large |fc|, the real part is 
proportional to k, with the constant of proportionality 
being equal to ±1 (within our numerical precision). This 
behavior is expected because at large frequency and mo- 
mentum, the vacuum of the boundary theory is effectively 
Lorentz- invariant. So, the dispersion relation of the exci- 
tations should effectively be relativistic at large frequency 
and momentum where the effect of charge density is neg- 
ligible. 

To have a better understanding of the spectrum, it is 
important to know how the residues (at the poles) be- 
have as a function of k. In order to numerically com- 
pute the residues of G-(w,k) at ui = cu^(k), we first de- 




FIG. 8: |Res G_ (u t , fc)| as a function of k for the leading negative- 
frequency pole in Figure 6(b) which is closest to the real axis and 
located to the left of the negative imaginary axis. We set d = 3, 
p = 5, q = 1 and m = 0. 



veloped series expansions for a_(z;cj,fc) and /3_(z;u;,fc) 
both near the horizon around z = 1 — e and near the 
boundary around z = 5. We then numerically integrated 
the (decoupled) differential equation for a-(z;uj,k) and 
/3_(z;w, k) from z — l^e to z — S and matched the nu- 
merically integrated solutions and their (first) derivatives 
with their boundary series expansions at z = 6. In so do- 
ing, we were able to compute the residues of G_(cj,fc) 
at cj = uj*{k) for a fixed k, denoted by Res G_ (u>* , k). 
Repeating the same steps for different fc's, one can nu- 
merically obtain the dependence of the residues on k. We 
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FIG. 9: A close-up of of the real (left plot) and imaginary (right 
plot) parts of the dispersion relation of the leading pole shown in 
Figure 6(b) which is closest to the real axis and located to the left 
of the negative imaginary axis. 



computed the fc-dependence of the residues of G_(u;, k) 
for the leading negative-frequency pole (closest to the 
real axis). Shown in Figure 8 is the absolute value of 
Res G_ (w* ,k) as a function of k for this pole. Also, the 
plots in Figure 9 show a close up of the real and imaginary 
parts of the dispersion relation of this pole. Comparing 
the plot of the residue to Figure 6(a), it is now appar- 
ent why the gap forms: as the leading negative-frequency 
quasi-normal mode approaches u> — 0, its residue dies off 
quickly. Since no other mode has appreciable spectral 
weight, this accounts for the suppression of the spectral 
weight near w = 0. For larger values of p, the residue 
falls off more rapidly, and consequently the gap widens. 



V. FINITE TEMPERATURE 

So far, our analysis has been at zero temperature. 
However, there are important aspects of Mott insulators 
that transpire at finite temperature. In particular, there 
are Mott insulators [20] exhibiting a transition to a con- 
ducting state as the temperature is increased. The classic 
example of this is VO2. Below = 340AT, VO2 becomes 
insulating with a gap of A = 0.6 eV. This ratio of the 
gap to the critical temperature A/T* is approximately 20. 
This behavior should be contrasted with systems such as 
superconductors in which U(l) symmetry is broken and 
A/T c w 1 — 2. That A/T* well exceeds unity is one of 
the unresolved puzzles with VO2. It points to strong 
correlations being the source of the gap rather than the 
breaking of some spontaneous symmetry as in the case of 
superconductivity. Optical conductivity studies [30] re- 
veal that spectral weight as far away as 6 eV contributes 
to the formation of the Drude peak at zero frequency 
once the Mott gap closes. Such UV-IR mixing is a ubiq- 
uitous feature of Mott systems. While we have argued 
that our holographic setup can capture the high-low en- 
ergy spectral weight transfer, we have not yet addressed 
the finite temperature aspects of the Mott problem. 

The boundary theory we are investigating here can eas- 
ily be studied at finite temperature by considering the 
RN-AdS4 background away from extremality, namely for 
< Q < y/3. Using the same procedures outlined above, 
we obtained the spectral function and studied the den- 
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FIG. 10: A close-up of the density plots of Im G— (uj, k) for p = 6 
and T/n ~ 5.15 X 10" 3 (left) and T/fi ~ 3.98 X 10" 2 (right). A 
gap is still seen in the plot on the left while it is closed in the plot 
on the right. 



sity of states as a function of temperature. As Figure 
10 reveals, the Mott gap observed here does in fact close 
as the temperature increases. Further, the transition is 
sharp. To estimate the ratio of the zero-temperature gap 
to the temperature at which the gap closes, T*, we take 
a close-up of the density of states and study its evolu- 
tion as a function of temperature, see Figure 11. Indeed 
for p = 6 (or p = 7), we find that A/T* ~ 10. Though 
smaller than A/T* in vanadium oxide, it does illustrate 
that dynamically generated gap we have found here does 
possess non-trivial temperature dynamics. 

For the record, we show in Figure 12 the quasi-normal 
frequencies of a_ (bottom plot) as well as ImG_(w,fc) 
(top plot) for k = 2, p = 5 and T/fi = 0.16. Notice 
that the branch cut at zero temperature has dissolved at 
finite temperature into a series of isolated poles on the 
negative imaginary axis 3 . Also, as it can easily be seen 
from the plots in Figure 12, the (real part of the) two 
quasi-normal frequencies of a_ which are closest to the 
real axis match quite well with the location of the peak 
on the left hand side and the bump on the right hand side 
in ImG_(w,fe = 2). The other quasi-normal frequencies 
represent poles in Im G_(u;,fc = 2) which are relatively 
wide and have small residue. Although not shown, to 
the extent that we have checked, all the poles stay on 
the lower half w-plane as k is varied. 



VI. DISCUSSION 

We have studied extensively the dichotomous behavior 
of the boundary theory fermion correlators in the pres- 



3 Note that in order to generate the bottom plot in Figure 12, the 
expression in (34) should be appropriately modified to reflect 
the fact that the system is at finite temperature. This is tied to 
the fact that at finite temperature f(r) has a single zero at the 
horizon. 
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FIG. 11: A close-up of the density of states A(lu) at p = 6 for 
T/fi ~ 0.44 (dotted), 0.16 (dashed) and 5.15 X 10~ 3 (solid). 
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FIG. 12: The plots in (a) and (b) show, for k = 2 and T/fi = 0.16, 
Im G— (u>, fc) as a function of u) and the quasi-normal frequencies of 
a- , respectively. Here, d = 3, p = 5, q = 1, m = 0, and M = 250. 



ence of a bulk Pauli coupling in our holographic set up. 
For the boundary theory dual to the extremal RN-AdS4 
background, we showed that as we vary p from large nega- 
tive values up to a small positive value of p = 1 /y/6 (while 
keeping m = and q = 1 fixed) , the behavior of the exci- 
tations change from Fermi liquid like (for p < —0.53), 
though not in the precise Landau sense in which the 
width of the excitations is quadratic in frequency, to a 
marginal Fermi liquid at p — —0.53 and on to a non- 
Fermi liquid for —0.53 < p < l/v6- In the context of 
the earlier work [8] in which such behavior was observed 




FIG. 13: A cartoon of the zero temperature "phase diagram" in 
the m, q, p parameter space. Different regions of the phase diagram 
correspond to each of the principal structures in the cupratc phase 
diagram (compare to Figure 1). 



by changing the scaling dimension, as well as the charge, 
of the boundary theory fermion operator, the Pauli cou- 
pling offers a more direct connection with Mott physics. 
Our argument here is based on the fact that for large 
positive values of p a Mott gap arises (as evidenced by a 
vanishing of the quasiparticle residue shown in Figure 8) 
in the spectrum of the boundary theory fermion operator 
without the apparent breaking of a continuous symme- 
try. This is Mott physics. We have seen these features 
by holding m and q fixed while varying p. It is clear 
that the basic properties that we have seen will persist 
throughout a domain in the m,q,p parameter space. In 
Figure 13, we suggest that some locus through the param- 
eter space can be identified with doping in the cuprate 
phase diagram: each of the principle features in the nor- 
mal state of the cuprates is present. It is interesting to 
compare this heuristic phase diagram with that of the 
cuprates in which there is a continuous evolution from 
a Mott insulator in the undoped state to a Fermi liq- 
uid in the overdoped regime. In between these extremes 
lie non-superconducting non-Fermi liquid states charac- 
terized by a pseudogap (a suppression of the density of 
states without any long-range superconductivity) and a 
strange metal in which the resistivity is a linear func- 
tion of temperature. Our work suggests that the Pauli 
coupling mimics the role of the electron filling. 

It is natural to investigate how the introduction of a su- 
perconducting condensate would complement the physics 
that we have discussed here. A suitable charged back- 
ground at zero-temperature was studied in [31], follow- 
ing [32, 33]. For a range of parameters (namely, for 
— 2q^ < —3/2 where and q^ are the mass and 
the charge of the bulk scalar field, respectively), a bulk 
solution with a non-zero charged scalar condensate is pre- 
ferred over the Reissner-Nordstrom solution. (Indeed, for 
this range of parameters, the Reissner-Nordstrom solu- 
tion is unstable against turning on the scalar field in the 
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FIG. 14: Density plot of the boundary theory fermion spectral 
function for (a) p = and (b) p = 3. Here, q^ = 1.5, L = 1, and 
£t = 2\/3. The black lines depict the IR light-like region and the 
red curves represent the bound states. 
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FIG. 15: kp as a function of p. Here, q^ = 1.5, L = 1, and 



bulk.) This solution is asymptotically AdS4, and = 
has a near horizon geometry that is also AdS4 with a fi- 
nite speed of light, ctr, and a finite dynamical exponent. 
The horizon of the zero-temperature solution (which is 
the Poincare horizon of AdS4) is at r = and there is no 
residual entropy at zero temperature. Although we will 
discuss this elsewhere [34], preliminary studies of the ef- 
fect of the Pauli coupling on the boundary theory fermion 
correlators 4 indicate that the main contributor to the gap 
in the fermion spectral density in the superconducting 
phase is the so-called Majorana scalar-fermion coupling 
(denoted by rj 5 in [36]). As there arc limited tools avail- 
able for the study of "electrons" at strong coupling, the 
model proposed here could offer key insight into how su- 
perconductivity emerges from a background in which all 
energy scales are coupled. As a first step, we set 775 = 0, 
and briefly discuss here the effect of the Pauli coupling on 
fermion correlators in a boundary theory dual to the su- 
perconducting background of [31]. Suppose = 0, and 
q$ > v / 3/2, so that the near horizon geometry is AdS4 
(with a characteristic radius Lir) and assume there ex- 
ists a (2+l)-dimensional IR CFT dual to this AdS4 near- 
horizon geometry. For definitcness, wc set q^ — 2q, where 
q is the charge of the bulk fermion 5 . The Dirac equation 
for ip± (as well as the corresponding flow equations for 
£±) and the IR boundary conditions can easily be worked 
out. A crucial difference compared to the case of the RN- 
AdS4 is that here the dimension of the IR CFT operators 
dual to if)±{r —> 0) does not depend on p. Depending 
on the sign of s 2 = — w 2 /c 2 R + k 2 , the Green functions 



4 See [35-40] where the authors analyze fermion correlators (in 
the absence of the Pauli interaction) in some superconducting 
backgrounds. 

5 This condition is not required when r?5 = 0. We consider this 
condition so that our analysis here can be generalized to the case 
where 775 [34]. Also, note that the convention of charge in 
[31, 36] is different than our convention in previous sections by a 
factor of two, namely ghere = 2gthere- In this discussion, we use 
the convention of [31, 36] for q and q^. 



G±(oj,k) exhibit different behaviors. (Note that one still 
has G_(w, fc) = G + (uj, —k).) For s 2 > (IR space-like 
region), the horizon boundary conditions for £± are real 
and since the flow equations are also real, one concludes 
that the boundary theory fermion spectral density is zero 
in this region, except when there are bound states (by 
which we mean poles of Re G± (w, k) in the s 2 > region) 
of the Dirac equation. For s 2 < (IR time-like region), 
on the other hand, the horizon boundary conditions are 
complex resulting generically in a non- vanishing fermion 
spectral density. The IR space-like region is the primary 
feature of the finite ctr theory that distinguishes it from 
the RN-AdS4 theory. In particular, there is no analogue 
here of the log-oscillatory region (at u; = 0), and thus 
no analogous mechanism for the disappearance of zero- 
frequency poles as p is increased. Indeed as we describe 
below, as p is increased, the zero-frequency pole moves 
but persists. Thus in the presence of superconductivity 
where the U(l) symmetry is spontaneously broken, the 
gap in the fermion spectral density is determined by the 
Majorana scalar-fermion coupling 775. This implies that 
the parameter responsible for Mottness is distinct from 
those involved in superconductivity. 

Figure 14 shows a density plot of the fermion spectral 
function for p = (left) andp = 3 (right), where q$ = 1.5 
and /j, = 2-^/3- The density plot for p — 0, which has 
been previously obtained in [36], is also shown for the 
purpose of comparison with the density plot for a non- 
zero value of p such as p — 3. As the plot in Figure 
14(b) shows, turning on a non-zero value of p suppresses 
the spectral density of the incoherent excitations (those 
in the IR time-like region) and pushes them away from 
the boundary (the IR light-like region, depicted by solid 
black lines). But, since at u> = there is a Fermi peak 
(and indeed an infinite number of long-lived bound states 
for other values of uj in the IR space- like region) , turning 
on p does not result in the formation of a gap. In fact, 
for the above-mentioned parameters, we find that, for 
p = 3, kp 6.0. For p = 0, the Fermi peak, for the 
same parameters, is at kp ss 1.5 [36]. Turning onp moves 
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around the location of this Fermi peak as shown in Figure 
15. 

We note again that there are regions of parameter 
space (in particular, q$ and m^) where the preferred 
geometry is either Reissner-Nordstrom or the supercon- 
ducting geometry. Our results indicate that while the 
fermion gap in the superconducting geometry is con- 
trolled by the Majorana scalar-fermion coupling [31], 
a Mott gap can still form in the Reissner-Nordstrom 
regime. The decoupling of these two effects is promis- 
ing in the context of the cuprate phase diagram. 
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